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Summary. — The problem of l-dimensional ultra-relativistic scattering of 2 iden- 
tical charged particles in classical electrodynamics with retarded and advanced in- 
teractions is investigated. 

PACS 03.20 - classical mechanics of discrete systems. 
PACS 03.50 - classical field theory. 



D , 

(~| , Introduction 

One of the unsolved problems in theoretical physics is Hamiltonian and quantum 
K^ , description of discrete relativistic systems with retarded and advanced interactions. Such 

^ • systems result from exclusion of field freedom degrees in field theories: while expressing 

5t ! fields by sources in classical field theories, or performing integration by fields in generating 

functionals of quantum field theories (carried out at definite boundary conditions) finite- 
dimensional theories arise, describing the motion of sources. 

Field exclusion in classical electrodynamics leads to Wheeler-Feynman system [1] — 
a relativistic theory with finite number of freedom degrees, whose action is Poincare- 
invariant functional of world lines of charges. This theory is equivalent to classical 
electrodynamics under the following boundary condition: field emitted by a charge is 
completely absorbed by other charges, so that outside some sphere (of radius greater 
than the Universe one) the resulting field of charges vanishes. The given formulation of 
electrodynamics has the following properties: 

1. In Wheeler-Feynman electrodynamics a field created by a point charge is a half-sum 
of retarded and advanced potentials (dislike the standard approach, which uses only the 
retarded potentials). Consequently, the time reversion symmetry in this theory evinces 
not only in the action, but also in solutions of equations of motion. 

2. T-nonsymmetric effects like radiative damping arise in multibody problems and 
have thermodynamic nature [1]. 

3. As the classical theory has finite number of freedom degrees, one should expect 
that the corresponding quantum theory has no divergencies. 
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Canonical quantization of Wheeler- Feynman electrodynamics has not been performed, 
because Hamiltonian formulation of this theory was unknown. Moreover, Lagrange equa- 
tions of motion in Wheeler- Feynman electrodynamics are not ordinary differential equa- 
tions, but belong to a poorly investigated class of functional equations (differential equa- 
tions with deviating arguments [2]), that has no elaborate approaches to solution: neither 
general analytical methods, nor steady numeric techniques. Existence and uniqueness of 
solutions for this class of equations have not been studied thoroughly. That's why a gen- 
eral structure of solutions of classical electrodynamics in Wheeler-Feynman formulation 
is still unknown. 

Particular results obtained in this domain are: solution uniqueness theorem for large 
distances between charges [3] ; exact solution of the problem for circular motion of charges 
in the theory with retarded and advanced interactions [4] , and the theory with retarded 
interaction and radiative reaction [5]; exactly solvable modification of the problem [6], 
with one charge influenced by retarded interaction, and the other — by advanced one; nu- 
merical solution of 1-dimensional 2-body problem [7], extended to velocities v < 0.9545c. 

Here we will continue the solution of the last problem to greater velocities. Besides 
that we will find new solutions in the range oi v < 0.9545c (along with those in [7]). 
Hamiltonian formulation of the theory constructed in [8] is used for the problem solution. 

The paper has the following structure: the first section describes the methods under 
the use, the second one presents the obtained results. In three appendices boundary ef- 
fects and limiting regimes are examined, and also a list of additional questions interesting 
for further investigation is given. 



1. Methods 

0. We will consider one-dimensional scattering of two relativistic particles of equal charges 
and masses. The methods used allow to study scattering of particles of different masses, 
however, the case of identical particles is of the most interest due to additional symmetries 
the problem acquires. 

We use the following system of units: light velocity c = 1, classical radius of a particle 
e^/mc^ — 1 (all distances are measured in classical radii). 

1. A method described in [8] is used for the problem solution. The clue idea of the 
approach is a choice of special parametrization of particles world lines, in which the 
problem gets the following formulation (fig.l). Let's consider a system of particles x„, yn, 
moving in 2-dimensional space-time so, that they are always located in vertices of a 
polygonal line assembled of light rays. It's enough to set a motion of particles x, and 
define y trajectories as J/^(t) = x+(t), yni^) = x~_^i{t) (in light coordinates x^ = x^ ± 
x^). The motion is defined by a system of differential equations for x^{t) (simultaneous 
by r), which is equivalent to the original one with deviating arguments. Then additional 
conditions are imposed onto the particles trajectories, ensuring their sewing together 
to the whole paths of x and y particles. Therefore, this approach reduces the problem 
to a set of ordinary differential equations, that defines particles evolution, and a set of 
equations for initial conditions and evolution time of the form F{X) — 0, that provides 
trajectories sewing together. 
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Fig.l. Light stairway parametrization. 



It's possible to restrict calculations to a finite number of trajectories a;„, n — 1...N. 
This possibility follows from exponential increase of distance between particles with large 
n: \xn — yn\ ~ g^", 9 = (1 + v)/{l — v), V — the velocity of the particles in center-of- 
mass frame (CMF); and starting from some N, their interaction becomes negligible. As 
q ^ oo while u ^ 1, for large velocities few {N = 2,3) steps of the light stairway 
are sufficient to move the particles out of the interaction region. Thus, this method is 
especially convenient in ultra-relativistic case (for small velocities more equations should 
be considered). 

The motion is defined by a system of Hamiltonian equations for coordinates and 
momenta: 



(1) 

with Hamiltonian: 



dH + dH 

Pn = 
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The Hamiltonian is a polynomial function of r„ (linear for each r„ ). It's convenient to 
rewrite the equations completely in terms of (x, r), see [8]. 

The Hamiltonian is a Dirac's constraint [9], that briefly speaking means the following. 
As the Hamiltonian is conserved in evolution: 



H 



dH 

dXn 



dH 

dpn 



0, 



phase trajectories lie on the surface H ~ Const. In Wheeler-Feynman electrodynamics 
only those trajectories, that lie on the surface H = 0, correspond to a physically mean- 
ingful evolution. This constraint, causing the fact that not all momenta are independent, 
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arises due to parametrical invariance of the action and is typical for all rclativistic the- 
ories. E.g. in a theory of free relativistic particle analogous constraint has a form of 
a "mass shell condition" p^ — m^ = 0, identically satisfied on all trajectories, because 
Pf^ = mifj^/Vi^ by definition. 

The sewing conditions have the form 

(3a) <(T) = x::+i(0) n=l...iV-l, 

(36) pM(T)=K+i(0) n = 2...N-2, /i = +, - ; 

p+(T) = p+{0) - m/2 ■ {x+{0) -x+{Ojr\ 
p^_i(r) + m/2 • (x^(r) - x-^_,{T))-' = p-^io), 

for n = 2...N — 2 finite values of (x^,p^) coincide with initial values of these quantities 
for the next trajectory(^). The quantities {xf^,p!^)T, determined by solving differential 
equations (1), are functions of initial data {xl^,pl^)o and "time" T. So (3) is a system of 
nonlinear equations on initial data and T. 

It has been shown in [8] that conditions (3) alone do not provide smoothness of 
trajectories sewing. One more condition comes from the action minimum principle: 

(4) <(T) ^ <+,(0), < = (^||) - |f ' *" = 9n-9i ( I 

which is equivalent to smoothness of the y„ and yn+i trajectories sewing together. This 
condition should be imposed for some n — 1...N — 2, and together with (3) it will 
provide smoothness of sewing together all x„,j/„ trajectories (except of the boundary 
ones: 2:1^2, xn-i,n, see Appendix A). 

Thus, (4iV — 4) equations (3), (4) and H = are imposed on (4iV + 1) unknowns 
(x„ (0),p„ (0), T). Among the remaining 5 freedom degrees, 4 correspond to trivial trans- 
formations of solutions: 

• 2 translations of solutions: 

X —^ X + c ; 

• Lorentz transformations: 

x^ -^ x^c, x^ -^ x^ /c: p'^ -^ P^/c, p^ -^ pc; 

• reparametrization (shift along the trajectories): 

and one freedom degree corresponds to trajectories deformation in variation of relative 
velocity of the particles. 

2. The trivial freedom degrees can be eliminated, e.g. setting Xi{0) — 0, r^ (0) — Const 
— a value of this constant can be taken from Coulomb approximation (see below) and 
fixed in all further considerations. Applying translations and Lorentz transformation to 
a solution found, one can move it to another reference frame, e.g. CMF. 

(^) For n = 1,N ~ I two p„-conditions disappear, see [8]. Also, in the action transforma- 
tion (2)^(4) in [8], carried out for finite trajectories, off-integral terms appear, resulting to 
asymptotically vanishing corrections to the remaining two conditions. 
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In order to fix rcparametrization freedom degree, let's add one more equation: {x{ — 
xj~)(0) = {x^ — x^){T) in CMF, geometrical interpretation of which is clear from fig. 2. 
Amongst all of parametrization fixing constraints this one is especially convenient, be- 
cause it prevents a solution from "slipping off" along a trajectory and ensures that 
boundary points are distant from interaction region, both in the past, and in the future. 



Xn(T) 




Fig. 2. Fixing reparametrization free- 
dom degree. 



One more variable should be fixed in the problem, e.g. r^{0), it will control particles 
velocities. The expression for v (initial velocity of particles Xi,yi in CMF) has the form: 
V = 2{rf {0)ri {0))~^ — 1, see [8]. It's convenient to introduce a parameter /i = (1 — v)^^ 
with range [l,-|-cx3), and express r^{0) by it. 

Equation H = is linear for each r„ and can be solved explicitly (e.g. for ^2^(0)). 
Thus, with the changes made, we have the system of (4A^ — 4) equations F{^, X) = 
for (4A'^ — 4) unknowns X, depending on one parameter /i. 

Newton's method was applied to resolve this system. Coulomb (non-rclativistic) so- 
lution of the problem at v — 0.5 was used as a starting point, it has the form (in the 
designations introduced) : 



(5) 



^ (''') = TTT ch 2vT, X (r) — -—7 (2vr + ch2vTsh2vT), 



2;''(r„), r„=n-l-iV/2 n = l...N, 



1/2 

^1 n = l...N, ul 

x^ 



^n + l 



1/2 



n= 1...N -1, 



n = 2...iV, 



-Mat, 



Then after a few iterations precise solution was found (^). Then the /i parameter was 



(^) The condition |F| < 10 ^ has been used as a finish criterion; while far from solutions 
IFI ~ 10^..10^ 
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increased by A/i, and the solution found was used as a starting point for new value of /i. 
The step A/x was chosen adaptively: it was 

• decreased by 2 times, if convergency has been lost; 

• increased by 2 times, if solution has been found; 

so A/i has been automatically kept in an optimal region(^). 

3. With the help of these methods the solution was continued up to velocity v = 0.937. 
At this value of velocity Jacoby matrix of considering system degenerates (see fig. 3), and 
Newton's method becomes inapplicable. In fact, the solution extrapolation enables to 
"jump over" this point and continue the solution (0) behind it. However, vanishing of 
the Jacobian indicates a change of topology for the set of solutions, therefore the vicinity 
of this point should be studied by other methods. 







det J(v) 




4e+12 






i 






(-) 


1 h/ 


2e+12 


(0) 


--..^^^^ / 


/ / 

/ 



2e+12 




^^^"-^j^ 
















(0) \ 


4e+12 






- 



0.9 0.91 0.92 0.93 0.94 0.95 0.96 Fig.3. Jacobian dependence on velocity. 

Let's extract in the Jacoby matrix J^ — dFi/dXj of size K x K, K ~ AN — 4, a 
submatrix J of size {K — 1) x {K — 1), with determinant not- vanishing in the critical 
point. Place this non-degenerate block to bottom right matrix corner by renumcrating the 
variables and equations. Now, fixing (/i,Xi), it's possible to resolve [K — 1) equations 
Fi{^^Xi,Xj) = 0, i, j == 2...K, for {K — 1) unknowns Xj, using Newton's method. 
Non-degeneracy of J guarantees that the solutions found are isolated: for fixed (/i,Xi) 
in the vicinity of the solution Xj there are no other ones. Further, we get Xj{^,Xi) 
dependencies, carrying out step- by-step change of parameters (/i,Xi), similar to the 
described above. Substituting these dependencies to the remaining equation, we get a 
function of two variables f{fi,Xi) — Fi{fi,Xi,Xj{iJ,,Xi)), a behavior of which is shown 
on the fig. 4 (in terms of variables v, AXi = Xi — X^, where Xi — solution (0)). 

At w < 0.937 the equation f{AXi) = has a single solution. In the critical point a 
tangent to the graph at zero is horizontal (that's equivalent to dot J = 0). At u > 0.937 
the tangent slope reverses its sign and two additional solutions appear. Thus, in the 
critical point bifurcation of solutions takes place: one solution splits into three. 



(^) Other optimization methods are described in [10] 
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Fig. 4. Dependence f{v, AXi). In the vicinity of tlie critical point this function defines a 
surface, known in catastrophe theory as Cayley surface [11]. 
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Fig. 5. Trajectories shape. 



Positions of the additional solutions are determined using dichotomy. Taking v far 
enough from the critical point one can continue the problem solution using Newton's 
method. 
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2. The results 

1. Shapes of trajectories corresponding to found solutions are shown in fig. 5. While pass- 
ing critical point no irregularities of trajectories and graphs of velocities and accelerations 
are observed. 

For the main branch (0) trajectories are symmetric with respect to spatial and tem- 
poral reflections P : x^ ^ — x^ and T : x'^ ^ — .t^. 

For the additional branches (±) trajectories are not P- and T-symmetric, but trans- 
form to each other in these reflections. At the same time, these trajectories are PT- 
symmetric. 

This effect: violation of solutions symmetry for symmetrically stated problem is re- 
lated with non-uniqueness of solutions. Indeed, a symmetry of equations of motion (and 
initial data) under reflections means only that the reflected solution is also a solution 
(maybe another one), in other words, a set of solutions under reflections transits to it- 
self. Only if a solution at the given initial data is unique (e.g. as for ordinary differential 
equations), one can conclude that the solution coincides with the reflected one, i.e. is 
symmetric. Thus, violation of P, T-symmetries for solutions is speciflc to differential 
equations with deviating arguments, which can have multiple solutions for the same 
initial data. 

2. The dependencies of parameters drmtm on v are shown in fig. 6 (the definitions of 
parameters see in fig. 5: dm = mina;^ — maxy^ in CMF, tm — temporal distance between 
the points, where extrema are reached). 

The dependence dm(v) for the main branch of solution has minimum at velocity 
Vm — 0.956, corresponding to dm = 0.9075. This dependence (up to velocity value close 
to Vm) has been found in Andersen and von Baeyer work [7]. Overlaying the graph [7] 
to the one obtained here, we see their exact coincidence. 



Note. The following algorithm has been used in [7] for the problem solution. For a given path 
of particle x the forces acting on it from a particle y were calculated, assuming path of the y is 
a mirror reflected image of path of the x. Then the accelerations were integrated, corrections 
to X trajectories were found, and the described process was iterated. First of all, let's note that 
this algorithm considerably uses assumption about P-symmetry of the trajectories, therefore it 
is only capable to obtain symmetrical solutions corresponding to the main branch (0). Then 
in order to find a solution the described process was performed at a fixed dm value, which 
was diminished while moving to the range of greater velocities. The iterations converged till 
dm = 0.9077 (that corresponds to « = 0.9545) and diverged for less dm values. This divergency 
was caused by proximity of a minimum of dm{v) dependence — there are no solutions for 
dm < 0.9075. 



For V ^ 1 solution in the main branch tends to Hill's solution [12], for which the 
minimal distance between particles in CMF is dm = 1- The trajectories are polygonal 
lines with vertices in points (±0.5, ±0.5), see fig. 5. There are i5-like peaks on the acceler- 
ations graphs corresponding to these vertices (see top right graph in fig. 7), one of them 
corresponds to retarded, and another one — to advanced interactions. 

For (±) solutions dm, tm -^ for v ^ I. The trajectories in CMF tend to light rays 
emitted from the frame origin. 
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Fig. 8. On the top: (4-) solutions in various reference 
systems (dashed lines: v = 0.95, solid lines: v — > 1.) On 
the bottom: acceleration of particle x in the system 2. 




It's interesting to investigate limiting shapes of the trajectories in reference systems, 
where one of the particles has zero velocity at t — > ±oo. There are two such systems 
(see fig. 8), in one of them solution collapses to the frame origin, as in CMF. In the 
other system another regime is established. The trajectories have small fractures (breaks 
of smoothness), located in the points A, A'. These fractures are related with (5-like 
acceleration peaks, which correspond to the interactions of one type (in fig. 8 - retarded 
for X, and advanced for y). Interaction of another type has smooth dependence on time 
in ranges CA, A'C and defines trajectories shape in these intervals. It's possible to give 
analytical description of this shape, see Appendix B. 

3. Poincare-invariance of action [8] follows to appearance of Noether's motion integrals, 
i.e. to the conservation of 



(6) 



translations generators P — /_,Pi 

i 

and boosts generator M = > ^iPi ^ ^7p7 



In this approach momenta p^ are expressed via the rf^ values, which are uniquely defined 
by velocities x^ , and the {x^_^_i — x^ )~^ values, defined by distances between interacting 
particles, see (2), (5). The terms {x^_^^ — x^ )^^ correspond to interaction contribution to 
the integrals of motion, which is present in a standard approach [1] (emitted but still not 
absorbed field). The difference from the standard approach is that integrals of motion 
(6) comprise a sum over a sequence of points on the trajectories, linked by light stairway, 
which propagates both to the future and to the past of the system. 

Conservation of integrals (6) has been used for the control of solutions accuracy in 
numerical analysis, see [10] for details. 
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4. The solutions X(ii) found here were put in the Internet for the convenience of readers 
and further problem study. 

http : //viswiz . gmd . de/~nikit in/f w_data/nodel . html . 

Conclusion 

The analysis of 1-dimensional scattering of 2 charges of equal masses in Wheeler- 
Feynman electrodynamics has revealed splitting of solutions: the scattering is uniquely 
defined by asymptotic velocities of the charges for v < 0.937c; for v > 0.937c there 
are 3 different solutions. Mirror refiection symmetry is violated in the splitting: one of 
three solutions is P-symmetric, the other two are not, but transit to each other under 
P-reflection. For limits of solutions at w ^ c analytical expressions have been derived. 

Acknowledgments. The authors thank S.N.Sokolov and G.P.Pronko for useful dis- 
cussions. The work has been carried out at Institute for High Energy Physics (Protvino, 
Russia) and National Research Center for Information Technology (Bonn, Germany) and 
has been partially supported by RFBR 96-01-01273 and INTAS 96-0778 grants. 

Appendix A. 

Boundary effects 

1. Sewing (3), (4) leads to appearance of fractures on the boundary parts of trajectories. A 
physical reason of this effect is the following. 

The equations of motion [8] have been derived minimizing the action, defined for finite 
trajectories, see fig. 9. Due to this fact the interaction between particles y and x switches 
on/off instantaneously, potential of y particle field has breaks at the points X2,xn- These 
breaks correspond to 5-like Lorentz force F, which causes trajectories fractures at the points 
X2,xn- Breaks of velocity at X2,xn leads to breaks of second derivatives at y2,yN-i, and so 
on. Coefficient at (5-function in F, which determines amplitude of these breaks, is proportional 
to rj^^ and vanishes while »; — > oo (A*' ^ oo or w — > 1). In the examined range of velocities the 
effects caused by instantaneous switching of interaction appear to be small (a break of velocity 
in fracture points is Au < 5 • 10"** for v > 0.9, A« < 5 • lO"'^ for v > 0.95, see [10]). 



Xn(T) 




Fig.9. 
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2. Other schemes of interaction switching {e.g. continuation of trajectories by infinite straight 
lines) are able to remove the fractures. However, these techniques are less logical, as they do 
not follow from the action minimum principle for finite trajectories. Consideration of infinite 
trajectories in stairway parametrization is related with additional difficulties, see [8]. 

3. For A'' = 3 only the last 2 equations remain in (36); in (4) n — 1. The trajectories 2/1,2/2 are 
sewed together smoothly, fractures appear in sewing a;i,a;2,a;3 trajectories together. 

For N — 2 conditions (36) are absent, instead of (4) one should use a condition 



nUT)-uUO) = 



{x7,-x-){T) (4-4)(0)' 



Ml = -^, U2 = 



(4') 



— Xi and X2 trajectories are sewed together with a fracture, caused by action of <5-like Lorentz 
forces from bounds of j/i trajectory. 



Appendix B. 

Limiting regimes 

1. A problem of particles motion, inffuenced by interaction of one type {e.g. retarded for x and 
advanced for y) is exactly solvable [6]. Indeed, in light coordinates such a motion is described 
by differential equations, simultaneous by r = x^ , see fig. 10a. 




b) 





The equations have the form: 



{x^ 



y 



+^2' 



y = 



4i+y" 
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(a;^ 



7+^2' 



X =y =T, 



(in the equations (6) from [8] only terms corresponding to considered interaction are kept) . The 
solution is found by elementary methods: 



-{u)=A{j{u)+g{u)), y+{u)=A{f{u)-g{u)), x' {u) = y' {u) = Ah{u) , 



(B.l) 
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1-W2 • Y 

u 



1 1 -L T^ " 

i^ ^ 1 — u 



,11 1 + " 
-, J + f I'l^i ' 

I — w' I — u 



A 



9{u) 
(1 - v^'/^ 

4y3 



2v 



1 — v^ \ 1 — u 



-1<U<1. 



+ 



l~v' 



The trajectories are symmetric w.r.t. frame origin (are conserved under the transformation 



{x,yr 



<y, 



<^ u 



—u). Solution is given for CMF, v is asymptotic velocity of the 
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particles in CMF. General solution can be obtained from (B.l), applying arbitrary Poincare 
transformations: 



(B.2) 



{x,y)+ ^{x,y)+ -0 + 0+, {x,y)- ^ {x^yY /C + D' 



The solutions (B.l) have the following important feature; for i" — > 1 in zero velocity (laboratory) 
frame, depicted in fig. 10b, the trajectories have polygonal shape. Interaction localizes in a 
vicinity of fractures A, A' . Sewing this solution and P-reflected one together in points B, B' , a 
limiting Hill's solution for the problem with interaction of both types is obtained, see fig. 10c. 

2. Similarly one can obtain analytical description of limiting trajectories, presented in fig. 8. 
F*0 




b) 




c) 




Fig.ll. 



Let's consider extremal trajectories for interaction of the 1-st type (fig. 11a) and take into 
account contribution of interaction of the 2-nd type. The trajectories regions, where particles 
move at light velocity, give no contribution to the interaction of the 2-nd type, because this 
interaction is suppressed by the factor (1 — «)/(l + v) ^> Q for « — > 1 (see (6) in [8]). In the 
internal intervals CA, A'C' the interaction of the 2-nd type does not vanish, these trajectory 
regions have the shape, described by eq. (B.l), transformed by P-reflection (a;, y) -^ {y,x)^ and 
transformations (B.2). Outer intervals are straight and are sewed with internal ones smoothly 
in points C, C' , while in points A, A' there are breaks of slope, of the same amplitudes as those 
in fig. 11a. 

Thus, on internal intervals the solution has the form 



(B.3) 



x-=A{f-g)/C, y-=A{f + g)/C, x+=y+=AhC. 



(We have excluded translations D from transformations (B.2) and further consider solutions 
symmetric w.r.t. frame origin. A more thorough analysis shows that there are no solutions with 
D 7^ 0). Requirements of sewing with straight trajectory regions have the form: 



da;+ 
X {u) = l, X {u)=Q, --^iu) = l (point A in fig. lib) <^ 

ax 



(B.4) C = 1/A(v)h(u), f{v,u) - g{v,u) ^ 0, G - 



l-v' 



(1 - uu)2A(i.)2/i('u)2 



-1 = 0. 



The function x~{u) in the interval —1 < m < 1 monotonically increases from — oo to +oo, thus 
the second equation in (B.4) can be unambiguously resolved for u (with v fixed). Finding this 
solution numerically and substituting it to the third equation, we plot G{v) function graph, see 
fig.l2. 
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G(v) 



0.05 




-0.01 



0.9 0.92 



( + ') 



det J{v) 



-0.4 



-0.6 




0.6 



0.996 0.998 



Fig.13. 



The equation G{v) — has the following solutionsC*): 

1) « ^ 1 corresponds to Hill's solution. 

2) v^ 0.6087622455403174 {u = 0.867439322472601, C = 0.5948623884508324) corresponds to 
limit of solution (-I-). In order to transform this solution to laboratory reference frame fig. lie, 
let's define the coefficient 



^-iS-^"^) 



-1/2 



, where x (u) = —1 corresponds to point C on fig.llb. 



Solving the last equation, we find u = -0.8674393224726012, C = 3.2378655431429. Carrying 
out Lorentz transformations with coefficient C on the trajectories, we get limiting solution, 
depicted in top right of fig. 8. 

3) V = 0.9437848540619277 {u = 0.9998999343906364, C = 0.1703458543931005), the solution 
is presented in fig.13. In the reference system fig.llb this solution is close to Hill's one. Small 
(1%) discrepancies of these solutions in the interaction region are observed only in vicinities of 
points C, C' (also there is a small discrepancy of the solutions far from the interaction region, 
caused by difference of velocities Vx for these solutions). In the laboratory reference system 
fig. lie non-trivial limiting shape establishes for this solution as well (correspondent coefficient 
C = 34.5179309569148). For transition to reference systems, where Vx —f 1, coefficient C — > 



(^) Equations (B.4) arc satisfied with 10 ^^ accuracy on solution 2 and 10 ^"^ accuracy on 3. 
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is required. This, particularly, means that in CMF at ii — * 1 this solution collapses to the frame 
origin. 

Our hypothesis: for large v the solution (0) undergoes one more bifurcation. As a result 3 
solutions appear, one of them tends to Hill's solution at ii -^ 1, while the others to new found 
solution and it's P-adjoint. The form of det J(«) dependence (see fig. 13) conforms to existence 
of root near v = 0.9995. The presented techniques can not continue the solution to velocities 
V > 0.999, due to computational difficulties described in [10]. 

Appendix C. 

Proposal for further investigation 

The main result of our work consists in the discovery of a fact, that the system considered 
has the same number of freedom degrees as it's non-relativistic analog (a system of two particles 
on a line, interacting via potential forces). 

In both cases the shape of the trajectories on a plane {x,t) depends only on particles 
asymptotic velocities (in our problem is defined by it uniquely for v < 0.937 and three-way- 
ly for V > 0.937). Phase space of Hamiltonian theory in the studied problem is essentially 
4-dimensional, as well as in non-relativistic case. The phase space structure (both topological, 
and symplectic), however, for these cases considerably differs. The scheme for construction of 
reduced phase space for the considered problem looks like the following. 

The set of constraints (including H = and sewing requirements) extracts 5-dimensional 
surface in extended phase space {x,p)!^. This surface looks like S5 — Ti x Va x Ti, where 
Fi = X{ii) — the found branchy curve( ), Va — Poincare group, 7i — group of shifts along 
the trajectories. Reduction of symplectic form dp'^ A dx!^ from extended phase space onto this 
surface follows to symplectic form, degenerate along the 7i (solutions shifts are generated by 
Hamiltonian H, included into the full set of constraints). Factorization S5 by 71 (or imposition 
of gauges like r ~ xf) gives 4-dimensional phase space S4 = Fi x Va, where the reduced 
symplectic form defines Poisson brackets. 




^(^ 



Fig.l4. 



In further investigations of the given problem it's interesting to carry out: 

1. Calculation of Poisson brackets for independent variables and construction of canonical basis 
in the obtained Hamiltonian mechanics. 

2. A study of non-symmetrical problem mi 7^ m2- The fact that the branches X^'^'^'(fi) of 
found solution are sewed together in one point of phase space is related to the symmetries of the 
problem. For mi 7^ m,2 in a vicinity of critical point recombination and formation of disjoint 
branches of solution are possible (see fig. 14). 

3. Calculation of second variation of the action to investigate it's extrema. Necessity of such 
investigation consist in the following. Let's consider smooth function / : R" -^ R, which has 
fe > 1 extrema in R". In this case all extrema can not be minima, there should be extrema of 



( ) see http://viswiz .gmd.de/''nikitin/fw_data/nodel .html 
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other type among them (maxima or saddles) . One should expect that not all of found solutions 
of the problem correspond to mininmm of action, there should be extrema of other types as 
well. 

4. Investigation of system responses to external perturbations and analysis of causal properties of 
the system. As advanced interactions can lead to violation of causality principle, it's interesting 
to answer the following question: can an observer related to one of the trajectories detect the 
external influence on it before it takes place? In the case of causality violation one should 
estimate the magnitude of the effect. 
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